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The Gallavotti - Cohen Fluctuation Theorem (FT) impUes an infinite set of identities between 
correlation functions that can be seen as a generalization of Green Kubo formula to the nonlinear 
regime. As an application, we discuss a perturbative check of the FT relation through these iden- 
tities for a simple Anosov reversible system; we find that the lack of differentiability of the time 
reversal symmetry implies a violation of the Gallavotti - Cohen fluctuation relation. Finally, a brief 
comparison with Lebowitz - Spohn FT is reported. 



Despite many proposals have been advanced, a general theory of the steady state of dissipa- 
tive systems is still lacking. Nevertheless, it is a remarkable fact that under suitable hypothesis 
something can be said, and physical predictions can be made; for example, chaotic hypothesis 
(CH), [1], stating that for the purpose of studying macroscopic properties, a system exhibit- 
ing chaotic motions may be regarded as a transitive hyperbolic (that is, AnosovJ one, see [2], 
implies two remarkable results: the Gallavotti - Cohen fluctuation theorem, [3], [4], holding 
for reversible systems, and a formula describing the linear response of nonequilibrium systems, 
due to Ruelle, [5], [6]. These results have been and are still widely studied in the physical 
literature; see [7] (where a relation which inspired the FT was empirically discovered), [8], [9], 
[10], [11], [12] for instance. In this note we focus on the first on these two results; in particular, 
in section II we show that FT has some nice implication on nonlinear response (in [13] it has 
been already pointed out that FT implies the usual linear response theory, that is Green Kubo 
formula (GK) and Onsager reciprocity relations), while in section III we discuss a check of FT 
in a simple Anosov reversible system. Interestingly, a simple perturbative calculation shows 
that the lack of differentiability of the time reversal symmetry operator implies a violation of 
the Gallavotti - Cohen fluctuation relation. Finally, in section IV we show that in presence 
of a non differentiable time reversal symmetry an identity equivalent to the result proved by 
Lebowitz and Spohn in [14] is true. 

The connection between FT and nonlinear response consists in the fact that FT implies 
identities between correlation functions of physical observables in a nonequilibrium steady state; 
this has been pointed out first in [15], where the authors considered systems whose evolution 
was stochastic and ruled by a master equation, and independently in [16], in the context of 
deterministic systems satisfying CH. The results of [16] are presented in sections II and III of 
this paper. 

Before turning to our results, we spend a few words on some of the main features of Anosov 
systems; we refer the interested reader to [2] for a modern introduction to the subject. Consider 
a generic discrete dynamical system Xk = S'^xq (discreteness only causes technical problems, see 
[17], since S can be thought as the map arising from the Poincare section of a system evolving 
in continuous time), and assume that S is Anosov; informally, this means that given a point x 
the nearby points separate exponentially fast from x in the future and in the past, except when 
located on a surface Ws{x) [stable manifold) or Wu{x) [unstable manifold), respectively for the 
future and for the past. 

It is a well known result that Anosov systems admit an invariant measure the so - called 
SRB measure; in fact, given a sufficiently regular observable F[x) the following equality holds: 
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apart from a set of points of zero volume measure. It is a remarkable fact that the SRB measure 



2 



admit in principle an explicit representation, similar to the equilibrium Gibbs distribution. 

Notice that at equilibrium, that is when the system is stationary and non dissipative, the 
chaotic hypothicsis implies the ergodic hypothesis, in the sense that assuming CH the SRB 
measure is the Liouville one; but in general, when dissipation in present the SRB measure is 
singular with respect to the volume, that is it is concentrated on a zero volume set. 

To conclude, the SRB distribution verifies a large deviation theorem (see for example [17] 
for a proof of this statement for a special choice of F and in the more complex case of Anosov 

flows). In fact, consider the finite time average / = ^J2j=-^ F{S^x); then, it is possible to 

prove that there are values /i, /2 such that if [a, b] e (/i, /2) then Prob^^ (/ e [a, b]) ~ e^'^'''^^\ 
in the sense that 

hin ilogPro6^^(/e[a,6])= max Cf(/), (1-2) 

T^+oo r /e[a,fe] 

and Cf(/) is analytic and convex in (/i, /2)- 



II. FLUCTUATION THEOREM AND NONLINEAR RESPONSE 



At present time, no universally accepted definition of entropy for a dissipative system has been 
given. Nevertheless, the rate of entropy production is a well defined quantity, and is proportional 
to the work per unit time made by the thermostats on the system; the proportionality factor 
is the inverse of the temperature of the termostats (setting to 1 the Boltzmann constant). In 
particular, for a special class of thermostats, the gaussian ones, the entropy production rate 
corresponds to the phase space contraction, that is to minus the divergence of the equation 
of motion, see [18]; this fact can be taken as a general definiton of entropy production, if one 
assumes that the steady state of a large system is not aflfected by the details of the termostatting 
mechanism which ensures the existence of the steady state. 

In the case of a system evolving in discrete time, which is the case that we want to consider, 
the entropy production rate a is given by <t(x) = — log | det dS{x)\, where dS if the jacobian of 
the time evolution S. 

Now, let p be the adimensional average over a time r of cr{x), that is p = Tlij^-^ a{S^x), 

where a+ is the SRB expectation of (t{x), i.e. a+ = J ^i+{dx)a{x), and call C,{p) the large 
deviation functional of a, as defined in (1.2); assume that CH holds, and that the system is 
reversible, which means that there exists a differentiable isometry / such that I o S = o I, 
I o I = 1. Then, as proven by Gallavotti and Cohen, see [3] or [19] for a proof detailed from a 
formal viewpoint, the following result holds. 

Fluctuation theorem (FT): There is p* >1 such that for \p\ < p* 

a-p)=ap)-p(r+, (2.3) 



This result has an interesting corollary. Setting nr{q)dq = Probfi^{p S [q,q + dq]), define 
A(/?) as 

A(/3)= lim -log /e^(«-i)'"+'^7r^(g)dg ; (2.4) 
T— »+c!o r J 

clearly, A(/?) is related to C,{p) through a Legendre transform, that is 

C(p) = max {I3u+{jp - 1) - A(/3)) , (2.5) 

and moreover \{(3) admits the following expansion: 

~\~^ on on 

m = Y. E <'-(-M^*^-)-^(5*"-^-)>^ ^-E^"^' (2-6) 

n>2 ti,...,tn-i=— oo ' n>2 

where by < ... we denote the cumulant with respect to the SRB measure It is straight- 
forward to see that the fluctuation theorem implies an identity for the generating functional 
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A(/3), see [14], [15], [12] and [16] for instance: it follows that, as a consequence of the relation 
'!Tt{p) ~ e'^^'^+7r^(— p), valid under the hypothesis of FT, 

A(/3) = A(-l-/3)- <T+ (2/3+1) . (2.7) 

Notice that the generating functional of the cumulants is usually defined as (sec [14], [15], 
[12]) A(/3) = lim^^+oo \og{e~'^P'^+^), and with this definition the relation (2.7) is replaced 
by the more familiar A(/3) = A(l — /3); but the two definitions are equivalent, since A(/3) = 
— — A(— /3). Formula (2.7) translates immediatly in a relation for a+; in fact, (2.7) evaluated 

at /? = becomes: 

/I \n 

= A(-l)-a+^<7+ = ^C„i^. (2.8) 



n! 



n>2 

Assuming that the entropy production a{x) has the form (t{x) = J2i Gijf\x) +0(G^), where 

{Gi}, and {jf'\x)} are respectively the forcing parameters and the corresponding currents, it 
has been shown in [13] that the identity 

Y.GiG,dG.GMG=, -T.G^G^do.GMG=0 ' (2-9) 

which is nothing else that (2.8) at second order, is equivalent to the Green Kubo formula, 
stating that 

t — — QO 

where by < • >o we denote the expectation with respect to the invariant measure at zero forcing 
(which by CH is the Liouville measure). Hence, formula (2.8) can be seen as a generalization 
of GK formula to the nonlinear regime, being an identity for (t+ valid for G ^ 0. Moreover, by 
taking derivatives with respect to (3 in the r.h.s. and in the l.h.s. of (2.7) we find that 

^^^^(-1)^+^ r.>2, (2.11) 

fe>0 

which is a nontrivial identity for the cumulants valid at G 7^ 0. Finally, these identities can be 
considerably extended by using a generalized version of FT. Consider a generic observable O 
odd under time reversal symmetry, i.e. such that 0{lx) = —0{x), let w be its adimensional 

average over a finite time r, that is w = "^j^-i. 0{S^x) where 0+ is the SRB average of 
O, and call C,{PiW) the large deviation functional of the joint probability distribution iTr{p,w) 
of p and w; then, under the same hypothesis of FT the following result holds, as a special case 
of a much more general result in [20]. 

Generalized Fluctuation Theorem (GFT): There are w* > 1, p* > 1 such that for 

\p\ < p* and \w\ < w* 

({-p,-w)=C{p,w)-pa+ . (2.12) 



One can define the generating functional A(/3i,/32) of the mixed cumulants of O, cr in a way 
analogous to (2.4), 

A(/3i,/32)= lim ilog /e^(«-^)'"+^i+^(*-i)°+^=^7rx(g,t)dgdt ; (2.13) 

r— >+oo T J 

again, A(/3i,/32) is related to ({pjW) through a Legendre transform, and moreover it can be 
expressed as 

A(/3i,/32) = ^ J2 ^^Cn,m (2.14) 

k>2 m,n>0: m-\-n=k 
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where Cn,m = 9^j^dJ^^\{(3i, (32)\ ^^^q, that is Cn,m is given by a sum over times of mixed cumulants 

of a{S*^x),0{S*' x), 0<i<n, 0<j<m. Then, in full analogy with what has already been 
discussed, it is easy to show that GFT translates into an identity for the generating functional, 
namely 

A(/?i, /?2) = A(-l - -P2) - (2/3i + l)a+ - 2P2O+ , (2.15) 
which implies the following relations: 



n = S^L})! V Ck- (21+1), 21+1 , . 

^+ 2'=-! ^ (21 + iy.(k - (21 + i)y. ^ ' 

k>2 l=0,2l<k-l ^ / \ \ II 



2k 



Ci,n-i = ^-^! (^k+i,n-i, n>2,l<n; (2.17) 

k>0 

equation (2.16) is obtained setting f3i = (32 = —\ in (2.15), while (2.17) can be proved differ- 
entiating with respect to /3i, (i2 the r.h.s. and the l.h.s. of (2.15). It is interesting to see what 
happens to formula (2.16) in the linear regime. Assuming that a{x) = X), Gijl°\x) + 0{G'^) 
we can rewrite (2.16) as 

^+ = E (fc-l)! 2^-1 E fc!^fe^^0.'= + O(G') 

k>2,keven k>Z.kodd 



2*^i'i+ E (fc-1)! 2^-1 ^i>'=-i+ E 

k>4,keven k>3,kodd 

= ^^M+ E ^^(f^0,.-^Ci,.)+O(G^), (2.18) 

fe>3, feodd ■ ^ ^ 

and from (2.17) we find that for k odd Co,k - ^(^i,k = 0{G'^). Hence, 

00 

0+ = -Ci,i + 0(G2) = 2 E < ^iS'x)0{x) >+ +0(G2) , (2.19) 

t= — oo 

which gives: 

^ 00 

5g.O+|g=o = o E < -^f >o . (2.20) 



2 

t=— 00 



Formula (2.20) describes the linear response of a generic observable odd under time reversal 
symmetry; this result can be seen as a special case of the remarkable linear response formula 
obtained by Ruelle, [5], valid in a much more general context. 



III. CHECK OF FLUCTUATION THEOREM FOR A SIMPLE ANOSOV SYSTEM 

In this section we will perform a check of the fluctuation relation (2.3) for a simple Anosov 
model, the perturbed Arnold cat map, starting from the idcntites (2.8), (2.11); as wc are going 
to see, the lack of differentiability of the time reversal simmetry operator implies a violation of 
(2.3). Consider the following discrete evolution on the bidimensional torus 

x^. = S^XjQ mod 27r , (3.21) 

where 

SeX = Sx + elix), S=(^\ ly (/2(f)) ' ^^-^^^ 
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and f2{x) a-re trigonometric polynomials. The map 5* is the so called Arnod cat map, 

which is the simplest example of Anosov map: in fact, the eigenvalues A-j., A_ of S are such that 
A+ > 1, A_ < 1. Moreover the map S is reversible, that is there exist / such that loS = S~^oI, 
I o I = 1 , where 




(3.23) 



Notice that since S is conservative (det S = 1) (t+ = 0, which makes meaningless the fluctuation 
relation (2.3) (the adimensional quantity p is not defined); but one can derive the analogous of 
(2.3) for the dimensional quantity p' = a+p, and this relation becomes trivial if the evolution 
is conservative because in this case p' = 0. 

Consider now e 0. By the structural stability of Anosov systems, the evolution generated 
by Se is still Anosov provided e is chosen small enough. In fact, for s < Eq a eonjugation H^, 
defined by the identity Se o He = Hg o S, can be explicitly constructed through a convergent 
power series in e, and it turns out that He{x) is Holder continues in x, [2]. In an analogous way, 
the SRB measure can be explicitly constructed, and it follows that the expectations of Holder 
continuous functions exists and are analytic in e, [2]. In particular it follows that, generically, 



a+=<-\og\dS,\>+>0 , 



(3.24) 



which means that the system is dissipative and so that the invariant measiirc is singular with 
respect to the volume; hence, the check of the fluctuation relation (2.3) is nontrivial in this case. 
Notice that the proof of FT, [19], requires that the evolution is reversible, and in particular 
that the time reversal symmetry operator is differentiable; in our specific case the existence of 
He implies that 1^ — H^ o I o H~^ verifies le o Sg = S^^ o 1^, but due to the mild regularity 
properties of H^, Is{x) is likely to be not differentiable. Then one can ask what happens to 
the fluctuation relation; to understand this, we make the choice fi{x} ~ sin(xi) + sin(2a;i), 
f2{x) ~ 0. An explicit computation shows that the linear response is still valid, as expected, 
and that 



Cs = -12e3 + 0(£4) , 



(3.25) 



i.e. formula (2.11) with n = 3 is false at the lowest nontrivial order in perturbation theory, 
since it tells that C3 = C4/2 + 0{e^), and C4 = 0(£^)- This is enough to say that (2.3) is 
violated; in fact, from (2.5) it follows that, [16], 



Cip) 



C2 



0-+ 
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[C3 + 0(£')]+0((p-l)V), (3.26) 



that is 



-C{p) + C{-p) = p<^+ + p 



C2 Cz 4. 



+ P' 



0{p^e') , (3.27) 



and (3.25) implies that the difference —Q{p)+Q{—p) is not hnear inp (the cubic term is nonzero). 
This result also shows that, as expected, cannot be differentiable: a check that would be not 
so easy without using the FT. 



IV. CONCLUSIONS AND COMPARISON WITH [14] 

The Gallavotti - Cohen Fluctuation Theorem implies nontrivial identities between correla- 
tion functions valid at non zero forcing G, which reduce to the usual linear response in the limit 
G — > 0; we have shown that it is essential that the time reversal transformation be (continu- 
ously) differentiable. By checking the FT relation through these identities, in a simple Anosov 
system where the time reversal symmetry operator exists but is not differentiable, we have 
found that the relation expected from the possible validity of FT cannot hold (while the linear 
response is still valid, as expected). Notice, however, that in physical applications the time 
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reversal symmetry operation is usually regular; for instance, it can correspond to the inversion 
of velocities (or to more subtle permutations of coordinates, sec [21]). 

To conclude, it is interesting to note that in presence of a non differentiable time reversal 
symmetry operator / a different FT holds, the Lebowitz - Spohn one, see [14]; strictly speak- 
ing this theorem has been proved in the context of general Markov processes, but it can be 
understood also in the case of deterministic chaotic dynamics. This result applies in particular 
to systems that are (suitably small) perturbations of reversible Anosov ones, since, as it has 
been pointed out in the previouos section, by the structural stability of Anosov dynamics the 
time reversal symmetry is not destroyed by the perturbation (although it will be in general only 
Holder continuous in x). 

The fact that a fluctuation relation still holds is a natural consequence of the Gibbian nature 
of the invariant measure describing the steady state; this has been pointed out in [22], where a 
FT for the one dimensional Ising model in an external field was derived, and then systematized 
in [23], where large deviation rules for Gibbs states involving transformations different from 
time reversal were discussed. 

Following step by step the proof of Gallavotti - Cohen FT, [19], it is easy to see that the 

large deviation functional ({p) of the dimensionless quantity p = Si=-^ a{S^x) where 
a{x) = —\u{x)+\u{Ix) and A„(x) is the sum of the positive Lyapunov exponents corresponding 
to the local unstable manifold Wu{x) of S verifies 

C{-p)=Cip)-p^+, (4.28) 

(of course under the physical restriction on p to vary within the analyticity interval of C, see 
[19]) which reduces to (2.5) if I is differentiable, since in this case A„(/x) = —Xs{x). Formula 
(4.28) is equivalent to the Lebowitz - Spohn fluctuation theorem; in fact, a{x) is proportional 
to the logarithm of the ratio of the SRB probabilities of the trajectories x^r^x^^+i, ...,x:r_i 
and of its time reversed, which is precisely the "action functional" introduced in [14]. 
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